A common model for superluminal propagation in absorbing atomic media and lossy 

metastructures 
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We analyse the superluminal propagation of narrow-band pulses at resonances in dissipative me- 
dia. The output waveform is an attenuated, undistorted, time-advanced version of the input which 
can be interpreted as the result of the interference of two scaled replicas of the input having a posi- 
tive relative delay. The proposed analysis applies both to the propagation in a passive bulk medium 
at an electronic resonance and in a dielectric waveguide coupled to a lossy micro-ring resonator. 
This latter case provides a clear insight of the underlying physical phenomena. 



Superluminal propagation is among the most striking 
phenomena associated to the propagation of electromag- 
netic waves. It has been known for long that it is possible 
to attain group velocities in excess of c when propaga- 
tion takes place in anomalous dispersive media A 
large number of systems have been studied, both theo- 
retically and experimentally, which hold group velocities 
larger than c. Garret and McCumber [2| considered a 
smooth pulse propagating through a dispersive absorbing 
medium and showed that pulse peak moves at the group 
velocity classically defined, even if greater than c or neg- 
ative. An experimental confirmation was given in [3]. 
Superluminal group velocity was studied in transparent 
media with inverted population by Chiao [4. Steinburg 
et al. investigated photon tunneling time [5[ and super- 
luminal propagation with negative group velocity in a 
medium with a gain doublet [6j. Bolda et al. considered 
propagation with negative group velocity due to a nearby 
gain line. [7| Later, Dogariu et al. demonstrated superlu- 
minal pulse propagation by transparent linear anomalous 
dispersion created through the use of two close Raman 
gain peaks [8| , and Schweinsberg et al. observed superlu- 
minal (as well as slow) pulse propagation in an erbium- 
doped fiber (EDF) 0. More recently, Jiang et al. have 
given an experimental confirmation that superluminality 
in an active Raman medium can be brought about by a 
single single frequency pump field. [l(| 

Although superluminal behavior might seem unphys- 
ical, it is always found to be the result of some sort of 
artifact, so that the group velocity being greater than c 
or negative at certain frequencies does never imply an 
affront to relativity. In , a frequency- varying concept 
of group velocity is chosen that maintains the meaning of 
the function du> / d[He(k)] even for pulses that experience 
considerable distortion through propagation; an experi- 
mental demonstration using this description is presented 
in 12j. From a mathematical point of view, absortion 



or gain resonances are responsible for the shaping of the 
linear dispersive properties of the medium which permit 
to produce spectral regions with associated superluminal 
or negative group delays in bulk dielectric media. On 
the other hand, superluminal propagation in metastruc- 
tures such as those based on coupled microring resonator 
(CMR) has also been predicted [13], so the question im- 
mediately arises of whether some connection exists, and 
to what extent, between both phenomena. As we shall 
presently see, the basic physics of the superluminal prob- 
lem can be very insightfully addressed through the analy- 
sis of a simple CMR structure having a lossy ring waveg- 
uide. We will show that the model of interference of the 
input pulse with scaled self-replicas not only explains the 
superluminal operation of the CMR and similar struc- 
tures, but can also be translated to the context of prop- 
agation through dielectric atomic (bulk) media near the 
absorption resonances. Indeed, superluminal linear prop- 
agation in such atomic media can be explained in terms of 
the interference of the input electromagnetic signal with 
the echoes produced by the resonant coupling with the 
medium via the linear polarization. We thus focus on the 
lossy CMR as a specially illustrative sample system. 

When only one mode of the electromagnetic (EM) field 
is involved, propagation problems can be addressed us- 
ing a scalar approach. Two representative examples are 
the propagation of plane waves in homogeneous media 
(the vectorial formalism in ll|, for example, is in prac- 
tical terms restricted to this case) and that of fields in 
single-mode waveguiding structures. The pulse propaga- 
tion problem with a plane wave is entirely equivalent to 
that of a waveguide mode, with the only proviso that the 
modal field profile across the waveguide is invariant for 
all frequencies in the pulse spectrum. This condition vir- 
tually applies in all cases of interest, even with ultrashort 
pulses containing very few cycles of the carrier. There- 
fore, in order to simplify the results, we use a general 
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approach using abstract scalar signals and linear systems 
for the discussion of the propagation problems. This does 
not restrict the validity of the results and shows that the 
same principles are appliable to any linear systems. 

We thus consider a generic superluminal system which 
may consist of either an atomic medium of finite length 
or a section of a suitable photonic system. We as- 
sume a time-localized propagating signal with a com- 
plex envelope (analytical signal) x(t). We call X(£i) the 
Fourier transform (FT) of x(t), with f2 = lo — loq and 
<x>o the carrier frequency. We will denote y(t) the com- 
plex envelope of the output signal, Y(Q,) its FT, and 
il(fi) = Y(p,)/X(Q,) the system response. Our study 
deals with superluminal propagation in resonant struc- 
tures, so we will assume that H(fl) has a resonance at 
£1 = 0. The "superluminal condition" is then defined as 
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(1) 



including the case r g < 0. For an atomic medium, tq is 
the group delay experienced during the propagation in a 
medium of the same thickness and the same properties as 
the one considered, but with the resonance at loq removed 
from its response. For a CMR such as that shown in Fig. 
3, To would correspond to the straight waveguide alone, 
without the ring waveguide. 

As seen below [c.f. expression ((4])], the superluminal 
effect in this case arises from a distortion caused when the 
propagation in the medium essentially produces a first 
replica of the input pulse with a delay To, and a second 
replica with a delay To + At. It is necessary that the 
input pulsewidth be greater than At so that the second 
replica can annihilate part of the first without significant 
distorison. As a result, the output pulse is apparently 
time-advanced with respect to the first replica. Calling 
h(t) =FT~ 1 [H(9)], i.e. the system impulse response, 
this simply means that 



h(t) = aS(t - t ) - bS(t - t - At), 



(2) 



where a and b are real positive constants. Taking the FT 
of Eq. ([2]) yields H(fl) — aexp(— iUro) — 6exp(— ifl[ro — 
At]). If the spectrum of x(t) is much narrower than the 
resonance peak, the phase is given by 



arg_ff(S!) ~ — SItq + arctan 



M7At 



(3) 



From expression (|3|), the relative group delay relative, 
defined as the delay increase wih respect to that of the 
medium or metastructure without the resonance, is given 
by 
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FIG. 1: Amplitude frequency response and relative group de- 
lay for the ring geometry described in the text. 



For superluminal propagation, it is necessary that b < 
a. Eq. Q also shows that, as a — > 6, the magnitude 
of the delay r g increases and the amplitude, given by 
|If(f2)| ~ a — b, decreases. 

We analyse a 1-ring "scissor" (side-coupled integrated 
spaced-sequence of resonators) as shown in Fig. 3. With 
an ideal lossless ring waveguide, it acts as an all-pass op- 
tical filter. However, the presence of radiation losses due 
to the curvature results in the appearance of a resonance 
in the transmission curve. Such losses play the same role 
as the absorption losses in an atomic medium, which, in 
the last analysis, are due to the electromagnetic energy 
radiated away by spontaneous emission. The response of 
the structure is given by [13j 
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where H(£i) relates the complex envelopes of the output 
and input modal electric fields, < 9 < 1 is the real 
transmission coefficient of the directional coupler (as- 
sumed lossless), I is the length of the ring waveguide, 
< a < 1 is the attenuation factor due to the radiation 
losses in the curved sections, and k — (w/c)n is the real 
propagation constant of the waveguide, with n being the 
modal index and lo = loq + Q. The carrier frequency is 
chosen at a resonance of the structure: loqI — m2it, with 
m an integer, thus kl = mlix + tlnl/c.We then see that, 
if a w 9 << 1 (but a < 9), expression (JSJ) simplifies to 
yield, by inverse transformation, the result l[2|): 



h(t) ~ 6S(t) -aS(t - At), 



(6) 



with At — Ih/c Thus, 9 = a and a = b. [Note that, in 
this case, To = as the straight, unloaded waveguide is 
immaterial in the model of Eq. 

We study the propagation of optical Gaussian pulses 
using the Transmission-Line Model method FJ| to solve 
the time-domain Maxwell equations in a two-dimensional 
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FIG. 2: Comparison of the normalized reference (A) and out- 
put (B) waveforms in the ring geometry for (a) 1.5, (b) 3, (c) 
6 and (d) 12 psec 1/e input Gaussian pulses. At ~ 1.6 psec. 



geometry. We compare two structures: A waveguide sec- 
tion and the same waveguide loaded with a micro-ring 
resonator built with two straight segments of length 15.25 
/jm joining two semicircles of 18 /im radius. This design 
permits a better control of the coupling from the waveg- 
uide to the ring. For all the (straight and curved) waveg- 
uide sections, the core refractive index is n\ = 3.361, the 
cladding refractive index is ni = 3.168 and the width of 
the core is 0.6 /zm. The separation of the guides in the 
coupling region is 0.25 /im. 

According to our model, the ideal conditions for the 
observation of superluminal propagation in the micror- 
ing resonator structure are defined by a large coupling to 
the ring so 9 is small and, simultaneously, large enough 
losses so a is comparable to (but smaller than) 9. We 
first inject a very short 5 fsec pulse in order to compute 
the transfer function. We compare the signal at a given 
distance from the ring output with the reference signal 
at the same plane propagating in the unloaded waveg- 
uide; these are considered the output and input signals, 
respectively. Figure [1] shows the amplitude response and 
the net group delay, computed from the phase response. 
The loss mechanism is supplied by the radiation in the 
curved sections which increases with frequency. The net 
group delay is negative, corresponding to superluminal 
propagation, up to a frequency limit when the loss in- 
crease sets 9 > a. From that point onwards, subluminal 
propagation at the resonances is found. As we approach 
the critical frequency, the magnitude of the net group de- 
lay becomes larger, in agreement with Eq. |T]), and the 
resonances become sharper. 

Figure [2] compares the normalized reference (A) and 
output (B) signal waveforms for four values of pulse du- 
ration: 1.5, 3, 6 and 12 psec. The carrier frequency of 
210.55 THz is tuned close to the center of the leftmost 
resonance in Figure [T] The round-trip time in the ring 



FIG. 3: Field amplitude distribution \E(x,z,t — to)\ at t 
when the leading tail of a 6 psec Gaussian pulse is entering 
the ring. 



is approximately 1.6 psec. We observe how, as the pulse 
spectrum narrows in relation with the resonance band- 
width, the propagation distortion disappears. The rela- 
tion between the output/input peak signal levels ranges 
from —16 dB in case (a) to —27 dB in case (d). Figure 
[3] shows the field amplitude \E(x,y,t — t )\ distribution 
at time to when the leading tail of the 6 psec pulse is en- 
tering the ring. For sufficiently long pulses, the shape of 
this distribution remains constant until the pulse finally 
leaves the structure, given a quasi-stationary picture. 

We now turn to the atomic media, where superlumi- 
nality occurs for narrow-band pulses near a resonance. 
The dielectric susceptibiliy can be approximated by 



X-Xo ~ Xr- 



(7) 



where xo (real and constant) is the background contri- 
bution to x from all the resonances above the one con- 
sidered, and A = 2f2/7 is the detuning factor normal- 
ized to the resonance width 7/2. We will now see that 
the concept of retarded interference still applies. The 
key point is that the role of the ring waveguide in the 
CMR structure is played by the electric polarization in 
the atomic medium. The atomic polarizability acts as 
an absorber and retarded reemitter of the propagating 
photons. Under certain conditions to be derived next, 
such function yields superluminal behavior analogous to 
that of the electromagnetic field when fed back into the 
straight waveguide through the coupler in the CMR. 

Starting from the wave equation for the electric field 
in the slab, E(z,t), denoting E(z,t) the complex field 
envelope, and using the slowly-varying envelope approx- 
imation (SVEA), one obtains 



dE(z,t) 



no 
c 



2a 

1 



dE(z,t) 
di 



-aE(z,t), (8) 



where tiq — (1+Xo) 1 ^ 2 is the background index of the ma- 
terial, a = Xr&/ (2jt.qc) is the (strong) amplitude attenua- 
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tion factor, with u> the carrier frequency, tuned to that of 
the atomic resonance. Eq. ([8]) is derived by considering 
the frequency dependence of the resonant polarization up 
to first order in lo — u), which is valid for sufficiently slow 
(spectrally narrow) pulses. From Eq. © it follows that 
the resonant group velocity is given by 
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2a 

7 
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(9) 



where the negative term modifies the background, 
nondispersive value c/no, and may lead to negative group 
velocities. There is a certain resemblance between this 
situation and the description of the lossless barrier tun- 
neling given by Winful [15J , where the apparent superlu- 
minal tunneling of narrowband pulses through a barrier 
is shown to be caused by the modulation of a standing 
wave. In our present dissipative, quasistatic problem, we 
see that the peak pulse does indeed not propagate super- 
luminality, but arises from an interference phenomenom. 

Now, using the transformations r = t— zuq j c and £ = 
z, an equation is obtained for the amplitude E'(Q,t) = 
E(C,,t + Quq/c) in the frame moving at the background 
group velocity. Further, except for the exponential atten- 
uation of the transmited pulse, the spatial and temporal 
variations are very slow: The global picture of the evo- 
lution of the optical field in the medium is analogous to 
that given in Figure for the CMR case. So, by writ- 
ing E'((,t) — exp(— a()f(C, t) the following equation is 
obtained for the slowly varying amplitude /(C, r): 



dj_ 

dC 



2adJ_ 

7 dr 



(10) 



In order to bring the result calculations into the for- 
malism of Eqs. ©-(HI), we approximate Eq. (fTO")) by a 
difference equation, which can be done because /(£, r) is 
a slow function, with the proviso that the slab thickness 
L is sufficiently small. Writing df((,r)/d( ~ [f(L,r) — 
/(0, r)]/L and df(C, r)/dr ~ [/(0, r) - /(0, r- At)}/ At, 
where At — Luq/c, we finally come to the result 



E(L, t) ~ aE(0, t) - bE(0, t - At) 



with 



2aL 



b = 



2aL 
7 At 



(11) 



(12) 



Expression (|11[) is entirely equivalent to expression ©, 
with x(t) = E(0,t) and y(t) = E(L,t). We see that 
\b\ < \a\ and \b\ — ► \a\ as (7At)/(2<xL) 0, which is the 
same regime as that considered for the CMR case. Fig- 
ure 4. a shows the input and output normalized waveforms 



for a 58 psec lie width Gaussian pulse propagating in 9.5 
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FIG. 4: (a) Normalized input and output waverforms for the 
superluminal propagation of a Gaussian input 58 psec pulse 
in a 9.5 /im thick medium with parameters obtained from Ref. 
[3|. (b) Output waveforms as obtained from the anaylitical 
solution and the difference equation (|lip . 



/tm absorbing slab with the medium parameters obtained 
from 0] • The output waveforms obtained using the ana- 
lytical expression E(L, t) — E(0, t — L/v g ) exp(— aL) and 
the difference equation (lll|) are displayed in Figure 4.b 
illustrating excellent agreement. 
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